Prospects for an electron electric dipole moment search in metastable ThO and ThF 
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The observation of an electron electric dipole moment (eEDM) would have major ramifications 
for the standard model of physics. Polar molecules offer a near-ideal laboratory for such searches 
due to the large effective electric field (F e ff), on order of tens of GV/cm that can be easily oriented 
in the lab frame. We present an improved method for simply and accurately determining P e ff, in 
a heavy polar molecule, allowing for a quick determination of candidates for an eEDM experiment. 
We apply this method to ThO and ThF + , both of which possess metastable 3 A electronic states. 
The values of F e fi in ThO and ThF + are estimated to be 104 GV/cm and 90 GV/cm respectively, 
and are therefore two of the best known candidates for the eEDM search. 

PACS numbers: 11.30.Er 
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One of the more spectacular goals of in the rapidly 
growing field of ultra-cold molecules [l[ is the search for 
the electric dipole moment of the electron (eEDM). If 
found, the eEDM would be a touchstone against which 
ideas beyond the Standard Model(SM), especially Su- 
persymmetry, could be checked. On the other hand, if 
it is not found but experiments can push the eEDM's 
value below its current experimental limit of \d e \ < 
1.6 x 10~ 27 e-cm[2j], then key supersymmetric extensions 
to the SM become much less likely to be trueQ. Thus, 
table-top experiments at the lowest achievable energies 
are a direct compliment to collider experiments at the 
highest. 

The current experimental limit on the eEDM origi- 
nates from a high-resolution spectroscopic experiment 
on atomic thallium [2]. The ability to make this kind 
of precise measurement in atoms stems from key in- 
sights by SandarsQ and contributions by Flambaum and 
Khriplovich^, Q. These works show that the apparent 
electric field acting on an electron inside an atom can 
be far larger than the field that is directly applied by a 
macroscopic laboratory apparatus. A further insight by 
Sandars[7| pointed out that, for an electron inside a po- 
lar molecule, the effective electric field F c ff can be even 
larger. This field is enhanced by relativistic effects and 
scales with nuclear charge Z as Z 3 , therefore preferring 
heavy elements, the same as in the atomic case. Effec- 
tive electric fields as large as tens of GV/cm have been 
anticipated in certain heavy molecules, the largest so far 
being in HgF where F off ~ 99 GV/cmQ. 

In practice, F e ff is manifested in the molecule by its 
ability to distort an electron orbital of nominal s symme- 
try on the heavy atom into a combination of an s- and 
p-orbital. For this reason, several candidate molecules 
are composed of a heavy (n s 2 ) 1 S-st&te atom with fluo- 
rine which is adept at drawing electrons toward itself - 
therefore accounting for the needed distortion. Molecules 
such as BaF@, YbF^, HgF and PbF@ have all been 



considered. Departing from this trend, DeMille et al. 
proposed the metastable a(l) state of PbO[ll]]. More re- 
cently, an alternative experimental effort has proposed 
to trap molecular ions, counting on the long trap life- 
times to improve the signal [12j. In response, several ions 
were proposed, including HI+ [l|, PtH+, HfH+0, and 

HfF+lUGJ]. 

The molecular ion experiment has stressed the desir- 
ability of molecules that are easily polarized in small 
laboratory electric fields. For ion traps in particular, 
this attribute is a necessity because strong electric fields 
would ruin the trap's characteristics. Molecules with this 
property that also have a reasonable value of F c ff, are 
of (scrd<5) 3 A electronic symmetry. In such a molecule, 
the overall A symmetry guarantees a small ^-doubling, 
hence easy polarizability, while the electron in the a- 
orbital may still contribute a large F ff[l4|. In addi- 
tion, the magnetic g-factor of a 3 Ai molecule is very 
small - a tiny fraction of a Bohr magneton because 
g.j ~ (.9lA — <7<?S). This is zero, apart from diamagnetic 
and radiative corrections 16], providing less sensitivity to 
magnetic field noise. 
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The ACME collaboration [17] has seized on the poten- 
tial utility of 3 A molecules in designing a new molecular 
beam eEDM experiment. This experiment will use ThO 
molecules, which possess a metastable 3 A state not far 
above the ground state. This molecule is also attractive 
inasmuch as several of its isotopomers have no nuclear 
spin, eliminating the complexities of hyperfine structure. 
Similarly, the JILA EDM team is interested in the iso- 
electronic cation ThF + , which is expected to be similar 
to ThO in its electronic structure, for use in an ion trap 
experiment [1 81 ] . 

Until now, the most reliable estimates of F e g were de- 
termined using elaborate relativistic many-electron cal- 
culations. In this Letter we instead develop a competitive 
method based on nonrelativistic molecular structure cal- 
culations perturbed by the Hamiltonian arising from the 
eEDM. The method, based on an initial approximation 
in Ref. [HI, is both fast an accurate. Indeed, we now re- 
produce, to within 25% the vales F e g for all species that 
have been treated by more elaborate relativistic theory. 
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This circumstance opens the door for researchers to per- 
form broad surveys of potential eEDM searches, and to 
evaluate the expected sensitivity of these experiments. 
As examples, we present in the Letter estimates of F ff 
for ThO and ThF + , finding them both to be excellent 
candidates. 

The simple idea behind the method is that the pres- 
ence of the eEDM causes a perturbation to the molecule's 
structure due to a relativistic effect [19] • In Dirac nota- 
tion the perturbation takes the form[5[ 
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Here a represents the electron's spin, assumed to point 
in the same direction as the eEDM; F is the local elec- 
tric field experienced by the electron inside the molecule, 
which is well approximated by the Coulomb field due to 
the large nucleus; and d e is the magnitude of the eEDM. 

The influence of this interaction on the molecular spec- 
trum is computed in perturbation theory. To do this, Ref. 
fl4| expanded the molecular orbital into atomic orbitals, 



|*moi) = e s \s)+e p \p) + 



(2) 



Here the coefficients e s and e p represent the contributions 
to l^moi) due to the s and p atomic orbitals of the heavy 
atom. Since the influence of the eEDM relies on relativ- 
ity, its contribution to the spectrum of the molecule is 
dominated by motion of the electron near the heavy nu- 
cleus. Thus the |s) orbital is essential. Moreover, since 
the eEDM Hamiltonian ([T]) has odd symmetry under par- 
ity, the matrix element (s\Hd\s) vanishes, and we must 
consider mixing with the heavy atom's \p) orbital. The 
other atomic orbitals that comprise l^moi) are far less 
significant to our purposes. 

The energy shift due to the perturbation is therefore 
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To evaluate AE, we estimate t StP using nonrelativistic 
molecular structure software. The relativistic effects oc- 
cur mostly in the matrix element (s\Hd\p). To evalu- 
ate this, we follow Ref. and use one-electron Dirac- 
Coulomb wave functions in place of the true orbitals |s) 
and \p) . This allows the integrals to be done analytically. 
Moreover, they are strikingly accurate provided that the 
quantum numbers n StP are replaced in the result by the 
effective quantum numbers v StPl to reflect the fact that 
these electrons actually exist around a structured ionic 
core, as opposed to a bare Coulomb potential. 
Using this approximation, the energy shift is(l4j 
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which identifies the item in square brackets as the effec- 
tive electric field. In this expression, a is the projection of 
the s-electron's spin onto the permanent dipole moment 
of the molecule; e is the electron's charge; and do is the 



Bohr radius, which means that the factor ej a\ identifies a 
characteristic electric field strength. The relativistic fac- 
tor, r ro i, accounts for the overlap of the Dirac Coulomb 
functions: 
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where 7 = y/(J + 1/2) 2 - {Za) 2 is a familiar dimension- 
less quantity, and Z is the atomic number of the heavy 
nucleus. By contrast, the effective charge Z e g is the 
charge of the heavy ion that the electron orbits. For 
example, in YbF, the bond is very nearly polar, and the 
valence electron contributed by the Yb + orbits about a 
Yb-core with Z c g — 2. Further, in Eq. (j4]), ho is a coeffi- 
cient (usually slightly less than unity) that accounts for 
the fact that |\P mo i) is typically a combination of several 
different orbital configurations; and xo is a further reduc- 
tion of order unity that accounts for spin-orbit mixing of 
l^moi) with other symmetries. Heavy diatomic molecules 
are often written in a basis where only the projection of 
the total angular momentum f2 is a good quantum num- 
ber. For example, the 3 Ai electronic level can have some 
admixture of a nearby 3 ni electronic level. 

The key to the success of the method is that the coeffi- 
cients e SjP can be extracted from nonrelativistic molecular 
structure software (specifically, we have used the molpro 
suite of codes [13] )• This is done via the overlap integrals 



.,p = (*mol| (s,p)), 



(6) 



where in Ref.(l4| the wave functions l^moi) and \(s,p)) 
were obtained from separate calculations, one for the 
molecule and one for the heavy atom. The atomic calcu- 
lation moreover identified orbitals with a particular prin- 
cipal quantum number n, associated with the atom's va- 
lence shell. In the context of distorting the heavy atom's 
s orbital, the appropriate p orbital must be the compo- 
nent p z , i.e., the projection with vanishing angular mo- 
mentum about the molecular axis. 

The true molecular orbital is a combination of many 
atomic orbitals. Therefore, the s-component is not nec- 
essarily a state purely composed of a principal quantum 
number n, though it is likely dominated by one of them. 
The presence of the light atom will distort the atomic or- 
bitals of the heavy atom and cause mixing of not only s 
and p but of n and n', where n' is another atomic energy 
level of the heavy atom. These other configurations will 
contribute to a science signal as well, albeit somewhat 
weakly. 

Yet, a basis expansion of \ij) mo i) can be done into any 
set of functions; for example, the Gaussian atomic basis 
set itself can be used. We can read off the contribution 
of each s and p heavy atom Gaussian to the molecular 
wavefunction. These are the functions that are varia- 
tionally optimized to produce the molecular wavefunc- 
tion and therefore utilize all the s and p character that is 
available to minimize the energy of the molecular state 
of interest. 
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Therefore, it is worthwhile (and more direct) to extract 
the coefficients e s p from the molecular orbital calcula- 
tion itself. A molecular orbital is comprised of a linear 
combination of Gaussian basis functions with the same 
projection of A. In the case of the c-molecular orbital, 
we can write out the expansion of the molecular orbital 



I*, 



E c £K>+E<K> 



E 



i i / 



E4 



pa 



(7) 



Jp 



In this expression each ket represents, not an atomic or- 
bital, but rather a Gaussian basis function. The coeffi- 
cients are therefore the direct numerical output of the ab 
initio MOLPRO calculation. The superscripts h (I) iden- 
tify functions centered on the heavy (light) atom. Each 
Gaussian function is moreover identified by whether it 
has s or p (or d, /, etc.) symmetry with respect to its 
atom. Viewed in this way, the net s-wave heavy atom 
character of |^ mo i) is given by the sum of all its s-wave 
contributions, meaning that 
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Here the denominator is needed because the individual 
Gaussian orbitals are not necessarily orthogonal to each 
other. A similar expression holds for e p . 

This method incorporates the influence of several re's, 
including those higher than the nominal valence orbital, 
which may be mixed in virtually. This advantage is also 
a drawback, however, because the lack of identification 
of a particular principal quantum number n leads to an 
uncertainty in identifying the effective quantum numbers 
v s , v in J5J). 

This same ambiguity was actually present in the pre- 
vious form of the calculation [l4j]. The atomic orbitals 
\(s,p)} from © merely represent the choice of a basis 
set. The basis could easily have been the neutral atom 
or the ion. Each choice of expansion would require a dif- 
ferent choice of v s ,p- Because the original Gaussian basis 
set is optimized for the neutral atom, we extract u s p from 
the neutral atom spectra. v s v are related to the atomic 
quantum defect n s . v through [2l| 



vi 



(9) 



where Ry is the Rydberg constant, and E n is the ioniza- 
tion energy of the heavy atom from the |n /) energy level. 
This method assumes that the spectrum for the atom fol- 
lows a Rydberg like series for higher-n levels. We then 
find the effective quantum number that arises from the 
structured core of the heavy atom. For example, E n for 
Yb in the ground configuration would be given as 

E n = E 6s = £([Xe]4/ 14 6s) - £([Xe]4/ 14 6s 2 )), (10) 



with [Xe] being the Xe core and the configuration 4/ 14 6s 
representing the first ionization level of Yb; the level with 
one of the 6s electrons ripped away. In a similar man- 
ner the quantum defect for the p level can be found by 
finding the difference from the 4/ 14 6s6p to the ionization 
state 4/ 6s. Similarly, for Pb we determine v s {vp) by 
considering ionization of the 6s 2 6p 2 ground state to the 
6s6p 2 (6s 2 6p) state of Pb + , respectively. 

Using the molecular orbital calculation directly to find 
e SjP also confounds somewhat the determination of Z e f{. 
If we had calculated atomic Yb, we would have Z c s — 1, 
whereas if we had used Yb + , we would have had Z c s — 2. 
Since we used neither, the value of Z e fi is not specified. 
However, the molecular calculation does identify the per- 
manent electric dipole moment of the molecule, as well 
as the bond length, via the distances r h (r l ) of the heavy 
(light) atom from the center of mass. Assuming charges 
q h and q l of the atoms, we can assert 



7 h h , I I 

rfmol = r q +r q 

Qmoi = q h + q l - 



(11) 

(12) 



For neutral (ionic) molecules, Q mo \ = 0(1), and we then 
solve for the individual charges q h and q l . Z e g on the 
heavy atom is then 



Zeff = q h + 1. 



(13) 



For ThO (ThF+) we find Z cS = 1.6(2.3), which we use 
in evaluating ((5]). 

In order to check the accuracy of this method, we 
applied it to several molecules studied previously and 
present the results in Table HI We used the same ab initio 
calculations as in Ref. 14], but extracted the values of 
e S)P as prescribed above. As is evident, the new method 
improves our older method, yielding results within 25% 
of previously published values, thereby making it use- 
ful for estimating experimental feasibility of candidate 
molecules. The largest deviation in the table is that of 
HfF + . This deviation can be attributed mainly to the 
peculiar ionization route of Hf, which goes from 5d 2 6s 2 
to 5d6s 2 , indicating that the d-electron interloper plays a 
large role in the low-lying Hf energy levels and our single- 
channel quantum defect approach is too naive. This 
deficiency would conceivably be cleared up by a multi- 
channel approach to the quantum defects, which we have 
not attempted here. 

These difficulties are not present for the Th atom. We 
therefore suspect that the predicted F ff in Table U prob- 
ably lie within 25% of the answer that the fully relativis- 
tic calculations will ultimately provide. In more detail, to 
construct the molecular wave function for Th we use the 
molpro software suite[2(|. The aug-cc-pVQZ basis of 
Dunning [2^] is used to describe the O and F atom's s — / 
orbitals. The ECP78MWB of the Stuttgart group[Il is 
used to describe the 78 electron core potential of Th and 
the aug-cc-pVQZ basis set to describe s-f orbitals [24]] . 
We performed the calculation with an occupied (active) 
space of {7,3,3,1} ({5,3,3,1}). We took several points in 
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TABLE I: Comparisons of published values of F e g to old re- 
sults from Ref. [1^ and new results in the present work. All 
values are given in GV/cm. 



Molecule 


-Published 


(Jldll4l 


INew 


BaF 


7A\9] 


5.1 


6.1 


YbF 


26 [10] 


43 


32 


HgF 


99[S] 


68 


95 


PbF 


-29 [8] 


-36.6 


-31 


a(l) PbO 


26.2[11] 


3.2 [221 


23 


HI+ 


0.34[13] 


0.57 


0.34 


HfF+ 


24[15] 


18 


30 


ThO 


N/A 


N/A 


104 


ThF+ 


N/A 


N/A 


90 



R separating Th from O or F between 2.75 and 4.5 ao, 
in order to determine the ground state bond length of 
r = 3.47 (3.73) a for ThO (ThF+) respectively. The 
value of r for ThO agrees well with the experimental 
value of r = 3.48 (25j. The final calculation included the 
present electronic state of interest ( 3 A) as well as X 
X U, 3 II, 3 E+ and *A for both ThO and ThF+. 

Based on this calculation we find that the effective 
electric field in ThO (ThF + )is approximately F c g w 
104 GV/cm (90 GV/cm). This similarity in F cff is per- 
haps not surprising since the two molecules are iso- 
electronic. These comparatively large values make the 
metastable 3 Ai states of ThO and ThF + competitive 
with even HgF as viable candidates for an eEDM search. 
There are two more parameters of experimental relevance 
that need to be estimated along with F c ff ; the ^-doublet 
splitting and the lifetime of the metastable state. 



An order of magnitude estimate of the f2-doublet split- 
ting can be estimated from the molecular parameters de- 
termined in the calculation and the use of perturbation 
theory 26]. We find that the splitting is on the order 
of 10" 6 cm" 1 in both ThO and ThF+ yielding a critical 
field of a few hundred mV/cm in the metastable states, 
which is the field required to polarize the molecule. 

The lifetime can be estimated by finding the amount 
of 1 IIi that is admixed into the 3 Ai state due to the 
spin-orbit interaction of the molecule. 1 FTi can decay to 
the ground 1 Eo state through electric dipole radiation, 
whereas the 3 Ai state would not decay via this route 
due to symmetry constraints. From our calculations, we 
find that the lifetime of the 3 Ai state is on the order of 
lms (lOOms)for ThO (ThF+). This is long enough to 
make a measurement of this state in a slowed molecular 
beam experiment [l7T] or in an ion trapfLsj]. 

In summary, we provide a modified and more efficient 
method of calculating F e g in polar diatomic molecules. 
Our determination of F e ff in several test cases lies within 
25% of fully relativistic calculations. We apply this 
method to the 3 Ai states of ThO and ThF+ and find that 
they are extremely good candidates for an eEDM search 
yielding effective fields of w 104 GV/cm and 90 GV/cm, 
respectively. In addition, the f2-doublet splitting in these 
molecules is very small, allowing for easy polarization. 
The lifetimes of the metastable states are long enough to 
make spectroscopic measurements [l7i [18} . 

This work was supported by the NSF. We thank the 
ACME collaboration and E. A. Cornell for useful discus- 
sions. 
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